Common energy scale for magnetism and superconductivity in underdoped cuprates: 
a yLiSR investigation of (Ca^Lai j.) (Bai.75_^Lao.25+a)Cu30j/ 
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We characterize the spontaneous magnetic field, and determine the associated temperature Tg, 
in the superconducting state of (Cai:Lai_a:)(Bai.75_i,Lao.25+i)Cu30j, using zero and longitudinal 
field ^SR measurements for various values of x and y. Our major findings are: (I) Tg and Tc are 
controlled by the same energy scale, (II) the phase separation between hole poor and hole rich regions 
is a microscopic one, and (III) spontaneous magnetic fields appear gradually with no moment size 
evolution. 



There is growing evidence that at low temperatures 
(T), cuprates phase-sepjarate into regions that are hole 
"poor" and hole "rich"El. While hole rich regions become 
superconducting below Tc, the behavior of hole poor re- 
gions at these temperatures is not quite clear. Some 
data support the existence of magnetic moments in these 
regions. In impure cases, like Zn or Li doped YBCO, 
the impurity creates-both the hole poor regionsa and the 
magnetic momentscfQ. In pure cases, such as LSCO, these 
magnetic moments are created spontaneously and un- 
dergo a spin glass like freezing at Tgd. However, there are 
still many open questions regarding these moments and 
the spontaneous magnetic fields associated with them. 
For example: Is there a true phase transition at Tg? 
What is the field profile and how is it different from, 
or similar to, a canonical spin glass? Is the field confined 
solely to the hole poor regions or does it penetrate the 
hole rich regions? Also, the interplay between magnetism 
and superconductivity is not clear. Is strong magnetic 
background beneficial or detrimental to superconductiv- 
ity? 

We address 
these questions by performing zero and longitudinal field 
muon spin relaxation experiments on a scries of polycrys- 
talline (Ca:jLai_:r)(Bai.75_:cLao.25+x)Cu30y (CLBLCO) 
samples. This superconductor belongs to the 1:2:3 fam- 
ily and has several properties that make it ideal for our 
purpose. It is tetragonal throughout its range of exis- 
tence < X ^ 0.5, SQ-jthere is no ordering of CuO chains. 
Simple valegce sumaS, more sophisticated bond-valancg 
calculationsQ, and thermoelectric power measurements^ 
show that the hole concentration is x independent. As 
shown in Fig. ^, by changing y, for a constant value of x, 
the full superconductivity curve, from the under-doped to 
the over-doped, can be obtained. Finally, for different Ca 
contents, parallel curves of Tc vs y are generated. There- 
fore, with CLBLCO one can move continuously, and with 
minimal structural changes, from a superconductor re- 
sembling YBCO to one similar to LSCO. 
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FIG. 1. Phase diagram of 

(Cai,Lai_a;)(Bai.75_2,Lao.25+a;)Cu30j,. The dashed lines in- 
dicate samples with equal Tg. Insets: Tc/T^'''' and Tg/T^^"" 
as a function of K{x)Ay where Ay — y — 7.15, and K{x) is 
chosen so that all Tc/T^^^ data sets collapse to a single curve. 

The preparation of the samples is described elsewherei. 
Oxygen content was determined using iodometric titra- 
tion. All the samples were characterized using X-ray 
diffraction and were found to be single phase. Tc is 
obtained from resistivity measurements. We also veri- 
fied usin|g-|TF-/xSR that CLBLCO respects the Uemura 
relational3 and that it is a bulk superconductor. 

The //SR experiments were done at two facilities. 
When a good determination of the base line was needed 
we used the ISIS pulsed muon facility, Rutherford Apple- 
ton Laboratory, UK. When high timing resolution was 
required we worked at Paul Scherrer Institute, Switzer- 
land (PSI). Most of the data were taken with a ^He cryo- 
stat. However, in order to study the internal field profile 
we had to avoid dynamical fluctuations by freezing the 
moments completely. For this purpose we used the "^He 
cryostat at ISIS with a base temperature of 350 mK. 
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FIG. 2. (a) ZF-fiSR spectra obtained max — 0.1, 
y — 7.012 sample. The solid lines are fit to the data using 
Eq. ^ the dashed line is a fit using the simulation as described 
in the text, (b) nSK spectra obtained in longitudinal fields 
from the x = 0.4, y = 6.984 sample at 350mK. (c) Polariza- 
tion curves generated by the simulation program as described 
in the text. 

Typical muon asymmetry [^(i)] depolarization curves, 
proportional to the muon polarization Pz(t), are shown 
in Fig. ^ (a) for different temperatures in the x = 0.1 
and y = 7.012 (Tc = 33. lA') sample. The change of the 
polarization shape with temperature indicates a freezing 
process, and the data can be divided into three temper- 
ature regions. In region (I), given by T > 8 K, the muon 
relaxes according to the well known Kubo-Toyabe (KT) 
function, typical of tiije case where only frozen nuclear 
moments are presentlld. In region (II), bounded by 8 K 
> T > 3 K, part of the polarization relaxes fast and the 
rest as in the first region. As the temperature is lowered 
the fast portion increases at the expense of the slow one. 
Moreover, the relaxation rate in the fast portion seems 
independent of temperature. Finally, at long time the 
asymmetry relaxes to zero. In region (III), where 3 K 
>; T, the asymmetry at long times no longer relaxes to 
zero, but instead recovers to a finite value. This value is 
~ 1/3 of the initial asymmetry Az{0). 

To demonstrate that the internal field is static at base 
temperature, the muon polarization was measured with 
an external field applied parallel to the initial muon spin- 
polarization. This geometry allows one to distinguish be- 



tween dynamic and static internal fields-. In the dynamic 
case the asymmetry is field independentllj. In contrast, in 
the static case the total field experienced by the muon is a 
vector sum of H and the internal fields, which are of order 

{B'^Y''^. For H > {B'^f'^ the total field is nearly paral- 
lel to the polarization. Therefore, in the static CclSG, clS H 
increases, the depolarization decreases, and the asymme- 
try recovers to its initial value. Because we are dealing 
with a superconductor this field sweep was done in field- 
cool conditions. Every time the field was changed the 
sample was warmed up above Tc and cooled down in a 
new field. The results are shown in Fig. ^(b). At an 
external field of 250 G, the total asymmetry is nearly re- 
covered. Considering the fact that the internal field is 
smaller than the external one due to the Meissner effect, 
this recovery indicates that the internal field is static and 
of the order of tens of Gauss. 

We divide the data analysis into two parts: high tem- 
peratures (region II), and base temperature. First we 
discuss region II. Here we focus on the determination of 
Tg. For that purpose we fit a combination of_a. fast re- 
laxing function and a KT function to the dataEj 



Az{t) ^ exp (-V)d \ + AnKT{t), 



(1) 



where Am denotes the amplitude of the magnetic part, 
A is the relaxation rate of the magnetic part, and An is 
the amplitude of the nuclear part. The relaxation rate of 
the KT part was determined at high temperatures and 
is assumed to be temperature independent. The sum 
Am + An is constrained to be equal to the total initial 
asymmetry at high temperatures. The relaxation rate A 
is common to all temperatures. The solid lines in Fig. ^ 
are the fits to the data using Eq. |l]. 
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FIG. 3. Tg vs. y. The horizontal solid lines are the equal 
Tg lines appearing in Fig ^. Inset: Magnetic amplitude as 
function of temperature for a x = 0.3 y = 6.994 sample. The 
arrow indicates Tg of that sample. 

The success of this fit indicates the simultaneous pres- 
ence of two phases in the sample; part of the muons probe 
the magnetic phase while others probe only nuclear mo- 
ments. As the temperature decreases Am, which is pre- 
sented in the insert of Fig. |[ grows at the expense of 
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An- At low temperatures Am saturates to the full muon 
assymettry. A similar temperature dependecne of Am is 
found in all samples. The origin of the magnetic phase 
is electronic moments that slow down and freeze in a 
random orientation. The fact that A is temperature in- 

1/2 

dependent means that in the magnetic phase 7^ (i?^) , 
where 7^ is the muon gyromagnetic ratio, is tempera- 
ture independent. In other words, as the temperature 
is lowered, more and more parts of the sample become 
magnetic, but the moments in these parts saturate upon 
freezing. 

Our criterion for Tg is the temperature at which Am 
is half of the total muon polarization as demonstrated 
by the vertical line in the insert of Fig. |[ The phase 
diagram that is shown in Fig. ^ represents Tg for various 
samples differing in Ca and O contents. This diagram is 
systematic and rather smooth suggesting good control of 
sample preparation. As expected, for constant x, higher 
doping gives lower Tg. 

We have singled out three groups of samples with 
a common Tg = 11, 8 and 5 K as shown in Fig. ^ 
by the horizontal solid lines. These samples are rep- 
resented in the phase diagram in Fig. |l| by the dotted 
lines. The phase diagram, containing both Tg and T^, 
is the first main finding of this work. It provides clear 
evidence of the important role of the magnetic interac- 
tions in high temperature superconductivity. In fact™ 
this phase diagram is consistent with recent theoriestil 
of hole pair boson motion in an antiferromagnetic back- 
ground. Those theories conclude that Tc oc Jus where 
rig is the superconducting carrierj-density, and J is the 
antiferromagnetic coupling energyfij. From the measure- 
ments at constant x we see that Tg oc Jf{ns) where 
/ is some decreasing function of n^. We assume that 
= ns[K(x)Ay] where Ay — y ~ 7.15 is chemical dop- 
ing measured from optimum, and if is a scaling pa- 
rameter which relates chemical to mobile-charge dop- 
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Pz{t) = 3 + 3 Pi\B\) cosij\B\t)B^dB 



(2) 



ing. Since T^ 



Jn^ (optimum), both TJT^^"" and 



Tg/T^^"" should be functions only of K{x)Ay. We find 
K{x) by making all Tc/T^^^ collapse onto one curve. 
This is demonstrated in the upper inset of Fig. |^. Using 
these values of K{x) we also plot TgjT^^^ as a function 
of K{x)Ay in the lower inset of Fig. Again all data sets 
collapse onto a single curve. This indicates that the same 
single energy scale J controls both the superconducting 
and magnetic transitions. 

We now turn to discuss the muon depolarization at 
base temperature. In this case all the muons experi- 
ence only a static magnetic field, as proven above. This 
allows one to reconstruct the internal field distribution 
out of the polarization curve. The polarization of a 
muon spin experiencing a unique field B is given by 
P^{t) = cos^{9) + sin^(6') cos(7|B|t), where is the angle 
between the field and the initial spin direction. When 
there is an isotropic distribution of fields, a 3D powder 
averaging leads to 



where p(|-B|) is the distribution of |B|. Therefore, the po- 
larization is given by the Fourier transform of p{\B\)B'^ 
and has a 1/3 base line. When the distribution of B is 
centered around zero field, p{\B\)B'^ is a function with a 
peak at {B) and a width A, and both these numbers are 
of the same order of magnitude [e.g. Fig ^(b)]. There- 
fore we expect the polarization to have a damped oscilla- 
tion and to recover to 1/3, a phenomenon known as the 
dip [e.g. the inset in Fig ^ (b)]. Gaussian, Lorcntzian 
and even exponential random field distributionlij, and, 
more importantly, all known canonical spin glasses, pro- 
duce polarization curves that have a dip before the 1/3 
recovery. Furthermore, a dipless polarization curve that 
staruates to 1/3 cannot be explanied using dynamical ar- 
guments. Therefore, the most outstanding feature of the 
muon polarization curve at base temperature is the fact 
that no dip is present, although there is a 1/3 tail. This 
behavior was found in all of pur samples with Tc > 
and also in Ca doped YBCOO and Li doped YBColla. 
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FIG. 4. (a) The internal field distribution extracted from 
the simulations for the case of correlation length ^ = 3 lattice 
constants, maximum moment size of 0.06/is and magnetic 
moment concentration p = 15%. Inset: The muon spin polar- 
ization for that distribution, (b) The same as above for the 
case of p=35%. 

The lack of the dip in Pz{t) can tell much about the 
internal field distribution. It means that (B) is much 
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smaller than A. In that case the oscillations will be 
over-damped and the polarization dipless! In Fig. ^ we 
show, in addition to the {B) ~ A case described above 
[panel (b)], a field distribution that peaks around zero 
[panel (a)]. Here {B) is smaller than A, and, indeed, the 
associated polarization in the inset is dipless. Thus in 
order to fit the base temperature polarization curve we 
should look for p(\B\)B'^ with most of its weight around 
zero field. This means that diverges like \/B^ at 

\B\ —f 0, namely, there is abnormally high number of low 
field sites. 

It also means that the phase separation is not a macro- 
scopic one. If it were, all muons in the field free part 
would probe only nuclear moment and their polarization 
curve should have a dip or at least its begining as in 
the high temperature data. The same would apply for 
the total polarization curve, in contrast to observation. 
Thus, the superconducting and magnetic regions are in- 
tercalated on a microscopic scale 20A){3. This is the 
second main finding of this work. 

The special internal field distribution, and the nature 
of the gradual freezing of the spins, can be explained by 
the intrinsic inhomogeneity of hole concentration. The 
part of the sample that is hole poor, and for that reason 
is "more" antiferromagnetic, will freeze, while the part 
which is hole rich will not freeze at all. The variation in 
the freezing temperature of different parts of the sample 
can be explained by the distribution of sizes aad hole 
concentration in these antiferromagnetic islandsc^. The 
large number of low field sites is a result of the fact that 
the magnetic field generated in the magnetic regions will 
penetrate into the hole rich regions but not completely. 

To improve our understanding of the muon polariza- 
tion, we performed simulations of a toy model aimed at 
reproducing the results described above. A 2D 100 x 100 
square lattice is filled with two kinds of moments, nu- 
clear and electronic. All the nuclear moments are of 
the same size, they are frozen and they point in ran- 
dom directions. Out of the electronic moments only a 
small fraction p is assumed to be frozen; they represent 
magnetic regions with uncompensated antiferromagnetic 
interactions. Since these regions may vary in size the 
moments representing them are random, up to a max- 
imum size. The frozen electronic moments induce spin 
polarization in the other electronic laoments surrounding 
them. Following the work of otherscll, we use decaying 
staggered spin susceptibility which we take to be expo- 
nential, namely. 



X (r) - (-1) 



exp(-r/^) 



(3) 



where r = UxOx + Uyuy represents the position of the 
neighbor Cu sites, a is the lattice vector, and ^ is the 
characteristic length scale. Because of this decay, at low 
frozen spin concentration, large parts of the lattice are 
practically field free (expect for nuclear moments). How- 
ever, the important point is that no clear distinction be- 
tween magnetic and field free (superconducting) regions 
exists. 



The muon polarization time evolution in this kind of 
field distribution is numerically simulated. The interac- 
tion between the muon and all the other moments is taken 
to be dipolar, and ^ is taken to be 3 lattice constantsHfl. 
The dashed line in Fig. | is a fit to the T = 350 mK 
data, which yield p = 15% and maximum moment size 
~ 0.06/iB . As can be seen, the line fits the data very 
well. However, as expected, the fit is sensitive to p£,^ 
only, namely the effective area of the magnetic islands, 
so longer ^ would have given smaller p. The field distri- 
butions and the polarization curve shown in Fig. ^ were 
actually generated using the simulation. In (a) the spin 
density is 15% while in (b) the density is 35%. 

In panel (c) of Fig. we show the spin polarization 
for different hole concentration, varying from 0% to 35% 
with the same ^ = 3. The resemblance between the simu- 
lation results and the muon polarization as a function of 
temperature in panel (a) leads us to our third conclusion 
that the freezing process is mostly a growth in the total 
area of the frozen AF islands. 

We are now in a position to address the questions pre- 
sented in the introduction. The appearance of sponta- 
neous magnetic field in CLBLCO is a gradual process. As 
the temperature is lowered microscopic regions of frozen 
moments appear in the samples, and their area increases 
but the moments do not. In the ground state the field 
profile is very different from that of a canonical spin glass 
or any other standard magnet. It could only be gener- 
ated by microscopic intercalation of an abnormal num- 
ber of zero field regions with magnetic regions without 
a clear distinction between the two. Finally, and most 
importantly, the phase diagram containing both Tc and 
Tg leads us to believe that these temperatures are deter- 
mined by the same energy scale given by J. 
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